Relaxation of Josephson qubits due to strong coupling to two-level systems 
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We investigate the energy relaxation (Ti ) process of a qubit coupled to a bath of dissipative two- 
level fluctuators (TLFs). We consider the fluctuators strongly coupled to the qubit both in the limit 
of spectrally sparse single TLFs as well as in the limit of spectrally dense TLFs. We conclude that 
the avoided level crossings, usually attributed to very strongly coupled single TLFs, could also be 
caused by many weakly coupled spectrally dense fluctuators. 
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I. INTRODUCTION 

With the current progress in fabrication, manipulation, 
and measurement of superconducting qubits it became 
crucial to understand the microscopic nature of the en- 
vironment, responsible for decoherence. Recent exper- 
iments showed considerable advances even in the pres- 
ence of the low- frequency, 1//, noise without eliminating 
its sources. The operation at the 'optimal point' in the 
quantronium [l| or the drastic reduction of unharmonic- 
ity in the transmon reduced dcphasing (increased T2) 
significantly. Especially in the transmon the relation 
T2 ~ 2Ti is now routinely reported, which implies that 
the decoherence is limited by the energy relaxation (Ti) 
process, i.e., the pure dephasing is suppressed. 

The microscopic nature of the dissipative environment, 
responsible for the energy relaxation, is still unknown. 
While the effects of the electromagnetic environment, in- 
cluding the Purcell effect, can be reliably estimated [1], 
the intrinsic sources of relaxation remain unidentified. 
There exist numerous indications that the charge and 
the critical-current noise are induced by collections of 
two-level fluctuators (TLFs) residing, possibly, in the 
tunnel junction or at the surface of the superconduc- 
tors [1, 0]. Several microscopic models of TLFs have 
been proposed [1, 0, B • It was suggested [1, [l^l that 
in charge qubits the relaxation is due to charge fluctu- 
ators which are simultaneously responsible for the 1/f 
charge noise On the other hand, the flux noise 

might also be due to a large number of paramagnetic 
impurities [l^. [isj. In addition, strong signatures of res- 
ident two-level systems were found, especially in Joseph- 
son phase qubits [lH . in which the Josephson junc- 
tions have relatively large areas. Recently the coherent 
dynamics of a qubit coupled strongly to a single TLF 
was explored with the idea to use the TLF as a naturally 
formed quantum memory [TgI . [i3| . 

Strong coupling to fluctuators as a source of decoher- 
ence of a qubit was the focus of research in the past, cf., 
e.g., Refs. llSl. [l9l . [20I . These works concentrated, how- 



ever, on single-level fluctuators, i.e., an electron jumping 
back and forth between a continuum and a localized level. 
Such a system maps onto an overdampcd dissipative two- 
level system 
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In contrast, here we study the effect 
of strong coupling to underdamped (coherent) two-level 
fluctuators. 

In this paper we analyze the properties of a qubit cou- 
pled to TLFs and, in particular, the relaxation process. 
On one hand, we study the dissipative dynamics of a 
qubit strongly coupled to a single fluctuator; these re- 
sults could be useful for the analysis of the experimental 
data, e.g., for phase qubits. On the other hand, we inves- 
tigate a qubit coupled to a collection of TLFs, and our 
findings allow us to speculate about some general fea- 
tures of the microscopic picture behind the phenomenon 
of qubit relaxation, which could account for the exper- 
imental observations. In the next section we discuss in 
more detail the analysis in this paper and possible fea- 
tures of the fluctuator bath. 



II. MOTIVATION AND DISCUSSION 

Below we investigate the dynamics of a qubit coupled 
to one or many TLFs. Our findings, together with certain 
observations concerning the experimental data, allow us 
to speculate about the properties of the collections of 
TLFs in real samples. In other words, we suggest a pos- 
sible structure of the fluctuator bath, which is consistent 
with these observations. 

In this paper, we first consider a single TLF coupled 
resonantly to the qubit, i.e., with its energy splitting close 
to that of the qubit. These results can be useful for the 
analysis of experiments, in which the qubit is coupled 
resonantly to a single TLF (while the other fluctuators 
are far away from resonance and do not contribute to 
the qubit 's dynamics) or strongly coupled to one TLF 
and weakly to the rest, which form a 'background'. In 
this regime the TLF strongly affects the qubit 's dynam- 
ics, and wc observe two effects: (i) coherent oscillations 
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with the excitation energy going back and forth between 
the qubit and the TLF; and (ii) the decay to the ground 
state due to the energy relaxation in either the TLF or 
the qubit. The oscillations themselves also show decay, 
dominated by dephasing processes. We describe the os- 
cillation and relaxation processes and determine the rel- 
evant time scales. 

Further, we discuss the dynamics of a qubit coupled 
to a collection of TLFs. Our motivation is based on 
the following observations from the analysis of the ex- 
perimental data: (a) strongly coupled TLFs {strongly 
coupled refers to a strong qubit-TLF coupling) were ob- 
served ex per imentally in phase qubits with large-area 
junctions [l4[l5l[T6j. In these qubits the Ti time shows 
rather regular behavior as a function of the energy split- 
ting of the qubit (in the regions between the avoided level 
crossings, which arise in resonance with the strongly cou- 
pled TLFs); (b) in smaller phase and flux qubits the Ti 
time often shows a seemingly random behavior as a func- 
tion of the energy splitting of the qub it \ n , [4|. |23| ; (c) the 
strong coupling observed in Refs. Il4l . llSl . Ila requires a 
microscopic explanation. For instance, a large dipole mo- 
ment of the TLF, ed, is needed to account for the data, 
where d is of the order of the width of the tunnel barrier 
and e is the electron charge; and (d) experiments [Tsj 
suggest a very high density of (spin) fluctuators on the 
surface of superconductors. 

Based on these observations we speculate about a pos- 
sible microscopic picture of the fluctuations, which could 
be consistent with these observations: First, one could 
expect in the analysis of the dependence of Ti on the 
level splitting that the contribution of each fluctuator is 
peaked near its level splitting (when it is resonant with 
the qubit and can absorb its energy efficiently). Further, 
one might assume that for a large collection of spectrally 
dense TLFs (that is with a dense distribution of the level 
splittings), the corresponding peaks overlap strongly, and 
the resulting Ti -energy curve is smooth (even though for 
a dense distribution the contributions of the TLFs are not 
necessarily independent). Indeed, this general picture is 
consistent with the data: in charge and flux qubits, with 
smaller- area junctions, the TLFs are not spectrally dense, 
and resonances with single TLFs can be resolved in the 
dependence of the relaxation rate on the level splitting. 
This may look as a seemingly random collection of peaks. 
In contrast, in phase qubits, with large-area junctions, 
there arc many TLFs (for instance, the TLFs could be 
located in the junctions so that their number would scale 
with the junction area); thus the spectral distribution 
of their level splittings is dense and almost continuous. 
This may produce a smooth Ti-vs. -energy curve. Fur- 
thermore, one can speculate about the structure of the 
fluctuator bath. Suggested scenarios of the microscopic 
nature of the fluctuators find it difficult to explain the 
existence of the strongly coupled TLFs, which were ob- 
served, for instance, in th e q ubit spectroscopy via the 
avoided level crossings [3, [IE [111 ■ 

In other words, in our picture each TLF is only coupled 



to the qubit, and they are essentially decoupled from each 
other. For each of them the coupling to the qubit is 
much weaker than observed in experiments [3, [TBI . [3] ■ 
For a dense uniform distribution of the TLFs splittings, 
usual relaxation of the qubit takes place. However, as we 
find below, if the level splittings of the TLFs accumulate 
close to some energy value (which may be a consequence 
of the microscopic nature of the TLFs), as far as the 
qubit's dynamics is concerned the situation is equivalent 
to a single strongly coupled TLF. Thus, in our picture, 
weakly coupled TLFs may conspire to emulate a strongly 
coupled TLF, visible, e.g., via qubit spectroscopy. Note, 
however, also the results of Ref. [Ij, |23j [Hi pointing 
towards single strongly coupled TLFs. 

To demonstrate this kind of behavior, we further study 
the regime, where two or more fluctuators are in reso- 
nance with the qubit. Our main observation in this case 
is that the fluctuators form a single effective TLF with 
stronger coupling to the qubit. 

For a collection of many TLFs with a low spectral den- 
sity, we estimate the statistical characteristics (by averag- 
ing over the possible spectral distributions) of the random 
relaxation rate of the qubit and estimate corrections to 
this statistics due to the resonances that involve multiple 
TLFs. 

Finally, we discuss collections of spectrally dense 
TLFs. In this case we identify two regimes. If the 
TLFs are distributed homogeneously in the spectrum, 
they form a continuum, to which the qubit relaxes, and 
the dynamics is described by a simple exponential de- 
cay. If, however, a sufficiently strong local fluctuation of 
the spectral density of TLFs occurs, the situation resem- 
bles again that with a single, strongly coupled fluctuator. 
This may explain the origin of the strongly coupled TLFs 
observed in the experiment. 



III. MODEL 

We consider the system described by the following 
Hamiltonian 

n n 

(1) 

The first term, the Hamiltonian of the qubit in its eigen- 
basis reads Hq = — ^ cr^, where Cq is the level splitting 
between the ground and the excited states, and cTj, is the 
Pauli matrix. Similarly, the Hamiltonian of the n-th TLF 
in its cigcnbasis reads Hf^n — — 2 ^f^^ ^2,n- 

We consider 

only transverse couplings with the strengths v±^„, de- 
scribed by the third term. We note that the qubit-TLF 
interaction (e.g., the charge-charge coupling) would typ- 
ically produce also other coupling terms in the qubits 's 
eigenbasis (longitudinal and mixed terms; cf. the discus- 
sion of the purely longitudinal coupling cx cr^r^ relevant 
for the dephasing by 1// noise, e.g., in Refs [H, \2T\). 
However, for our purposes (description of relaxation) the 
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transverse coupling is most relevant since it gives rise to 
spin-flip processes between the qubit and TLFs. In our 
model all TLFs interact with the qubit, but not with 
each other. This assumption is reasonable, since the 
TLFs are microscopic objects distributed over, e.g., the 
whole area of the Josephson junction. Thus, they typi- 
cally arc located far from each other, but interact with 
the large qubit. The term i^sath describes the coupling 
of each TLF and of the qubit to their respective baths. 
We model the environment of the qubit and of the TLFs 
as a set of baths characterized by the variables Xi and 
coupling constants Pi (specific examples are provided be- 
low). We write down and solve the Bloch-Redfield equa- 
tions [H, |2§| for the coupled system of qubit and TLFs. 

In the course of solving the Bloch-Redfield equations 
many rates (elements of the Rcdfield tensor) play a role. 
It is useful to reduce those rates, when possible, to the 
'fundamental' rates, i.e., those characterizing the decou- 
pled TLFs and the qubit. Each fluctuator is thus char- 
acterized by its own relaxation rate 7^'" and by the pure 
dcphasing rate 7^'", with the total dephasing rate given 



by 7r = (l/2)7r +7^'". 



We define these rates be- 



low and also discuss the generalization for the case of 
a non-Markovian environment. Also the qubit is char- 
acterized by its intrinsic (not related to fluctuators) re- 
laxation rate j1 and the pure dephasing rate 7^. Again 
72 = (l/2)7i +7^- In what follows the mentioned rates 
are treated as fundamental. All the other rates, emerging 
in the coupled system of the qubit and fluctuators, are 
denoted by capital letters F. 



IV. SINGLE TLF 



We flrst consider a system of a qubit and a single TLF. 



H 



1 1 

£9 cr^ 



Bath 



(2) 



We restrict ourselves to the regime ~ ey 3> f_L. Then 
the ground state |0) w \g t) and the highest energy level 
|3) ~ |e I) arc only slightly affected by the coupling. On 
the other hand the states |e t) and |g j) form an almost 
degenerate doublet. The coupling v± lifts the degeneracy 
to form the two eigenstates |1) = — cos | \g D+sin | |e t) 
and |2) = sin | \g |) + cos | |e t) (cf. Fig. [T|). Here we 
introduced the angle tan^ = v±^/5lj where 5lu = eq ~ ef 
is the detuning between the qubit and the TLF. The 
energy splitting between the levels |1) and |2) is given by 



A. Transverse TLF-bath coupling 

First, we consider the simplest case, in which only the 
TLF is coupled to a dissipativc bath and this coupling is 
transverse. The coupling operator in Eq. ^ takes the 



v±=0 



\ei)- 
|eT>,l5 i)- 
|5T>- 



i(0l9 i) + cos(^0|eT) 
»(f) l9J)+s'n(0 let) 



■|3) 
|2)i 

|1) J ^osc 
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FIG. 1: Level structure of the coupled qubit-TLF system in 
resonance Sui = 0. For v± = the middle levels form a 
degenerate doublet. The coupling lifts the degeneracy and 
splits the levels by the oscillation frequency uJosc- 



form 



-ffBath — 2^f.^ ' : 



(3) 



where the bath variable X/^j, is characterized by 
the (non-symmetrized) correlation function Cf^±(t) = 

(^Xf^j_{t)Xf_±{0)'^. In thermal equilibrium we have 

Cf^±{-Lj) = e^'^/'^Cf^A^{ijj). We assume here that T < 
i.e., that the temperature is effectively zero, so that 
we can neglect excitations. 

We solve the Bloch-Redfield equations [H, ^ for the 
coupled system using the secular approximation. As the 
initial condition we take the qubit in the excited state 
and the TLF in its thermal equilibrium state. Tracing 
out the TLF's degrees of freedom we find the dynamics 
of (O (Fig.©. 

<cr,(t)> 




t [1/ri] 



FIG. 2: ((Tz) as a function of time in units of the inverse 
TLF relaxation rate 7/ for the case of the qubit exactly in 
resonance with the TLF (dotted black line). One observes 
oscillations with frequency Uosc- The solid red curve gives the 
decay averaged over the oscillations, characterized by aav and 
Tav and the dashed blue curve shows the envelope described 
by aenv and Tenv Parameters in this plot are (in units of 7/ ): 
£9 = 6/ = 100, = 10 
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FIG. 3: Level-structure of the coupled qubit-fluctuator system 
in the simplest case where only the TLF couples to a heat 
bath. The rates and Fjo in Eq. (|8} lead from levels |1) 
and 1 2) respectively to the ground state |0). The excited state 
1 3) is not included in this illustration. 



For the expectation value {uz) we find the following 
expression 

-rose t ^ ' 



+aosc cos{uJosct)e 



We can introduce the single decay rate in two different 
ways, either effectively averaging over the oscillations or 
including all parts and describing the envelope curve (cf. 
Fig. [2). In the first case, averaging over the oscillations, 
we choose f{t) = a^^i e"^^'^ * + ai^2 e~^^-^* and obtain 
for the amplitude and the decay rate 



H,2 



2^' 



+ 



2±L 



(9) 



(10) 



This gives a quasi-Lorentzian line shape of Tav{Sio) with 
the width of the order of the coupling v± and the maxi- 
mum value at resonance of TaviSu; = 0) = i 7/. 



where {fy^)^ w — 1 is the zero-temperature equilibrium 
value. We can separate the right-hand side of Eq. ^ into 
damped oscillations, with decay rate Tosc, and a purely 
decaying part. The amplitude and the decay rate of the 
oscillating part are given by 



(5) 
(6) 



V 

1 / 
2^1 



2.0 r 




/ 0.5 - \'^o 



Here the rate 
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(7) 



is the relaxation rate of the fluctuator. We observe that 
the decay rate of the oscillations, Tqsc, is independent of 
the coupling strength v±_ and of the detuning 5uj. Note 
that the physics considered here is only relevant near 
the resonance e/ « Cg, and we assume that the spec- 
trum Cf ±{iS) is sufficiently smooth in this region, so that 
C/,±(e/)«C/.^(e,). 

For the purely decaying part we find 



ri.i 



0-1,1 



^ 10 



where vlr, and Vir, are the rates with which the states ll) 



ro 



2sin4^, 
2 

COS - 7( , 



2 cos^ I , 



ri,2 = r: 



.2 i 



7(, (8) 
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and |2) decay into the ground state |0) (cf. Fig. [3]). 

As we can see, the decay law for {(Tz{t)) — Wz)oo 
given by a sum of several exponents. It is sometimes 
useful, e.g., for comparison with experiments where no 
fitting to a specific decay law was performed, to define 
a single decay rate for the whole process. If a function 
f{t) decays from f{t = 0) = a to f{t 00) = 0, we can 
define the single decay rate F from f{t)dt = a/T. 




6uj/v± 



FIG. 4: Amplitudes a and rates F for the decay of the oscil- 
lating (dashed red) and purely decaying (solid blue) part of 
the qubit's (az) as a function of the detuning Suj between the 
qubit and fluctuator. The detuning is taken in units of the 
coupling v±. 



Figure |4] shows the amplitudes and rates characteriz- 
ing the decay of the oscillations (dashed red) and of the 
purely decaying part (solid blue) of the qubits (cTz)- 

To describe the envelope we choose f{t) = 
aui e-^i.^ * + a^,2 e'^^'^ * + a^sc e"^"" * (Fig. ^. This 



5 



gives 



where 7^ is the pure dephasing rate of the TLF: 



1 / 



2vl {vl 



(11) 



with amphtude aenv ~ aj,,i + + losc = 2. This again 
gives a quasi-Lorentzian peak with height TenviSto = 
0) = i 7^ and width similar to that of Eq. pU]) . 

The results presented above are valid in the regime 
when the coupling between the qubit and the TLF is 
stronger then the decay rates due to the interaction with 
the bath, v± ^ 7/. In the opposite limit the Golden 
rule results hold, with the relaxation rate ^ (cf. 
Refs. fil, id ISC' 



(15) 



Here 5*^^1(0;) = i [C/j| (w) + C/.|| (— tj)] is the sym- 
metrized correlator. 

Similarly, the rates due to the qubit 's longitudinal cou- 
pling to the bath are given by 





= cos2^7^. 




^ 12 




II ( — Wosc) , 


^ 21 




II i^osc) , 


lip 




0). 



(16) 



B. General coupling 

We now provide the results for the general case, when 
both the qubit and TLF are coupled to heat baths. The 
coupling in Eq. is given by 

+ i (/3,,|| a, ■ Xq^w + Pq^i_ ■ X,,i) . (12) 

It includes both transverse (_L) and longitudinal (||) cou- 
pling for both the qubit and the fluctuator. The temper- 
ature is still assumed to be well below the level splitting 
tq~ ef so that we can neglect excitation processes from 
the ground state. 

We specify now the main ingredients of the Bloch- 
Redfield tensor of the problems. As in Eq. ^ the re- 
laxation rates from the states |1) and |2) to the ground 
state due to the transverse coupling of the fluctuator are 
given by 



10 



2 « f 
cos^ I 7j 



F: 



sin^ I 7/ : 



(13) 



Similarly the transverse qubit coupling gives rise to new 
rates, 



Figure [H] gives an illustration of the processes involved in 
the formation of the Redfield tensor. 



4, r2i 



r 



20 



10 



■|o) 



FIG. 5: Illustration of the relevant transition processes in the 
general case of arbitrary coupling of the qubit and fluctuator 
to a heat bath. In addition to the transitions from the central 
levels |1), |2) to the ground state |0) with the rates Fio and 
r2o, we now also have transitions between the two central 
levels with the rates T12 and r2i. The excited state |3) is 
again omitted in the illustration. 



For {<Jz{i)) we again obtain the decay law, Eq. 
The amplitude of the oscillating part, aosc, is still given 
by Eq. ([5]). The decay rate of the oscillations is, however, 
modified: 



Tosc^ ^(ri+Fi2+F2i)+F, 



(17) 



^ 10 



sin^ i 7« 



^ 20 



cos^ 4 7^ , 



2 '1 ' ^ 20 ~ ^"'^ 2 



(14) 



The longitudinal coupling to the baths, oc CTz, Tz, gives 
two types of additional rates in the Redfield tensor, a 
pure dephasing rate, F^, and the transition rates between 
the states |1) and |2), 

F;^ = cos2^7^ 



J- 12 
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4 '^i'^^^ {-^osc) 



The rates without a superscript represent the sum of the 
respective rates for the qubit and TLF: 

ri-7/'+7?, F^ = F^+F9, 

ri2 = ri2 + ri2' r2i = r^i + ni- 

The purely decaying part is given by a slightly more 
complicated expression. Defining 



A = F10 + F12 
C = 



B — F20 + F21 



(A-B)"+4Fi2F2i , 
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we obtain 



- (l + cos^e) 

2 cos e (A - B) + (ri2 + V21) sin^ ^ 



2C 



{A + B±C) 



In the limit || = /3jn = [5q^\_ ~ 0, we reproduce the 
results of the previous section. 

The decay of the average is again characterized by 



0-1,1 + = 



fc2 



ri,2 



(18) 
(19) 



We work in the experimentally relevant limit lOqsc ^ T. 
Then we obtain 



r{2 = r{i = sin2^r{ 



p9 - 
^ 12 



r^, ^sin^cr* , 



where 



1 

4 
1 



r?, 



4 /^qJI '^9,11 (^O'sc) ■ 

The decay rate of the average then reads 



(20) 



2 {vl + (5tj2) 



(Fi + 4r^) + 47( 5a;2 



(21) 

where = T^ + r«. 

At resonance the resulting relaxation rate Tav is the 
mean of the decay rates of the qubit and TLF. Thus, if 
the TLF relaxes slower than the qubit (as was the case in 
Ref.[li,r,, decreases. Figure IH] shows the average decay 
rate Tav for the two cases with 7^ bigger (solid blue) 
and smaller (dotted red) than 7'. The double-peaked 
structure in the first case is due to the contribution r„ 
of the longitudinal coupling to the baths. Exactly in 
resonance and far away from resonance the effect of F^ 
vanishes, while for 5ui ~ v±_ it produces somewhat faster 
relaxation. 



Non-Markovian effects 

In the previous section we have shown [see Eq. p7p ] 
that pure dephasing affects only the decay rate of the 
osciUations of {uz) ■ The result, Eq. pT]) . is valid 
for a short-correlated (Markovian) environment. More 



specifically the relations F' 



PI, 



5„ II {oj = 0) and Tf 



cos^ ^ 7? with 7« 
cos2 ^ 7/ with 7^ 




FIG. 6; Tav as a function of the detuning 5uj (rates in units 
of the qubit 's relaxation rate 7', detuning in units of the 
coupling strength vx) for the general case when the TLFs 
relaxation rate 7/ is higher (solid blue, 7/ = 1.57j) / lower 
(dashed red, 7/ — 0.67') than the qubit 's relaxation rate 7J. 
The parameters in this plot are = r{ = 0.3, ux = 5 (in 
units of 7^*). 



\ P'j |[ iS/ji {ui = 0) are valid only in the Markovian case. 
The generalization of these results to the case of non- 
Markovian noise, e.g., 1// noise, is straightforward [l|]. 
Assuming that at low frequencies Sq^\\ = Aq \\/uj and 
S'y II = Ay_||/cij, we obtain 

-r' t ^ ' 



-aosc cos {ujosct) fi/f(t) e 



A, 



and 



where ln/i/y(t) ^^^2 ^ (^f3'^^^^Aq^\\ pj^^^^f^\\ 

Fggj, now includes only Markovian contributions. 

In resonance, Suj = 0, we have cos ^ = and therefore it 
seems at the first sight that the 1// noise does not cause 
any dephasing. Yet. as shown in Ref. [sil, in this case the 
quadratic coupling becomes relevant. The instantaneous 
splitting between the middle levels of the coupled qubit- 
TLF system, |1) and |2), is given by 



^^osc(A:q,||,X/J|) 



V± 



(23) 



This dependence produces a random phase between the 
states |1) and |2) and, as a result, additional decay of the 
oscillations of (ctz)- We refer the reader to Ref. |3l|for an 
analysis of the decay laws and times. Thus slow (1//) 
fluctuations make the decay of the coherent oscillations 
of {az) faster without considerably affecting the average 
relaxation rate Tav 

For strong 1// noise, thus, a situation arises in which 
the oscillations decay much faster than the rest of (uz). 
In experiments with insufficient resolution this may ap- 
pear as a fast initial decay from (cr^) = 1 to (cr^) 1— Oosc 
followed by a slower decay with the rate Tav ■ 



7 



C. Two two-level fluctuators 

As a first step towards the analysis of tlie effect of many 
TLFs, we examine now the case when two fluctuators are 
simultaneously at resonance with the qubit. This situa- 
tion in the weak coupling regime was considered, e.g., in 
Ref [l^l. In this regime the fluctuators act as indepen- 
dent channels of decoherence and thus the contributions 
from different TLF are additive. However, in the regime 
of strong coupling between the qubit and the fluctua- 
tors, which is the focus of this paper, we do not expect 
the decoherence effects of the two fluctuators to simply 
add up. This means that the resulting relaxation rate is 
not given by the sum of two singic-fluctuator rates. The 
Hamiltonian of the problem reads 



-Z <^f,nTz,n + -z(^x V ±^„Tx ,n+ J^Ba.th , 



n=l 



71=1 



(24) 

where i?Bath contains now the coupling of each of the 
fluctuators to its respective bath. In the regime of our in- 



terest, e„ 



e/,2 3> i,w^_2i the spectrum splits 



into four parts. The ground state is well approximated by 
\g tt)- Analogously, the highest excited state is close to 
\e II). The coupling 

is mainly relevant within two 
almost degenerate triplets. The first triplet is spanned 
by the states with one excitation: {|e ft) , \g |t) , \g ti)}- 
In the second triplet, spanned by {|.g H) , |e ti) , |e it)}, 
there arc two excitations. At low temperatures and for 
the initial state in which the qubit is excited and the fluc- 
tuators arc in their ground states, only the flrst triplet 
and the global ground state are relevant. Within the first 
triplet the Hamiltonian reads 




where the energy is counted from the ground state. 

First, we consider the two fluctuators exactly in res- 
onance with each other, e^.i = ef 2 — £/, and approxi- 
mately at resonance with the qubit: w ef. We per- 
form a rotation in the two-state subspace spanned by the 
states, where one of the TLF is excited, by applying the 
unitary transformation 



1 

[/ = I cos a sin a 
— sin a cos a 

Choosing the angle a = arccos 
the transformed Hamiltonian 



V 







>l.i + ^'i,2 





(25) 



= , we arrive at 



\ 



2e// 



Figure [7] gives an illustration of what happens. Af- 
ter the rotation [Eq. (P5)) ] the qubit is coupled to only 
one effective state |1), whereas it is completely decou- 
pled from the "dark" state |2). For symmetric coupling 
v±_i = v±_2 the states |1) and |2) are just symmetric and 
antisymmetric superpositions of \g jf) and \g "f j). Thus, 



leTT) Is it) Is Ti) 



k-TT) |i> 



|2) 



FIG. 7: Illustration of the performed transformation in the 
one-excitation subspace of the Hamiltonian for the case of 
equal couplings v±^i = dx,2 ~ v±. Applying the rotation ((25)) 
we arrive first in a situation where the state |2) (the dark 
state) is completely decoupled from the other two states. The 
renormalized coupling then splits the remaining two states, 
giving a situation analogous to the coupling to one TLF. 

the rotation demonstrates that the situation is equivalent 
to only one effective TLF coupled to the qubit with the 
coupling strength 



Analyzing the coupling of the effective TLF to the dissi- 
pative baths (of both fluctuators), we conclude that the 
effective TLF is characterized by the relaxation rate 



7f 



2 /.I , 2 /,2 



At this point we can apply all the results of Sec. lIVAl with 
v± replaced by the renormalized v±_ and the relaxation 
rate 7/ replaced by 7/. In particular, using the formulas 
(O and pU]) we introduce 



r(i)(<5c.) = r.„(5^)|^,^_,^^^,^^. 



(26) 



Here the superscript (2) stands for couphng to two flue- 
f 21 

tuators. The function Tav ^Suj) is peaked around Suj = 0. 
The height and the width of the peak depend on the 
relations between the coupling strengths and the relax- 
ation rates of the two fluctuators. In the limiting cases of 
'clear domination', e.g., for v±^i ^ v±,2 and 7/^ 3> , 
everything is determined by a single fluctuator. In the 
opposite limit of identical fluctuators, i.e., for v±^i = v±^2 
and 7/^ = 7/^, the height of the peak [Eq. p6|) ] is given 
by 7i ^ /2 = 7i^ /2, exactly as in the case of a single fluctu- 
ator. The width of the peak is, however, \/2 times larger 
since v±_ ~ V2vj_^i = ^/2v±_2- Clearly, the relaxation 
rate of the qubit is not given by a sum of two relaxation 
rates due to the two fluctuators. 

If the fluctuators are not exactly in resonance, this 
result still holds as long as their detuning Sef = e/,i — e/,2 
is smaller than the renormalized coupling v±. For much 
larger detuning, the rate is given by the sum of two single- 
TLF contributions. 
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D. Many degenerate fluctuators 

For a higher number of fluctuators in resonance, i.e., 
e/^„ — ^/j the argument presented above is still valid, 
and the resulting decoherence of the qubit's state is the 
same as for a single TLF with a renormalized coupling 
strength of 




and an effective TLF relaxation rate 
-/ ^ 2 

n 

It should be stressed that this equivalence holds only 
within the one-excitation subspace of the system. 

If the system with many fluctuators is excited more 
than once it is no longer equivalent to a system with one 
effective TLF. Multiple excitation could be achieved e.g., 
by following the procedure used in Ref. [Ill or that of 
Ref. repeatedly, i.e., exciting the qubit while out of 
resonance, transferring its state to the TLFs, exciting the 
qubit again and so forth. The simplest case is when all 
the fluctuators have equal couplings to the qubit v^^n = 
v^. The system's Hamiltonian reads then 

H = -]^f'qCz - e/Sz + v±axSx , (27) 

where Sa = (1/2) J^n '^a,n- 

For procedures of the type used in Ref. [3, [11], i.e., 
when only the qubit can be addressed, the TLFs will 
remain in the spin representation of Sa in which they 
were originally prepared. If the TLFs are all initially 
in their ground states, the accessible part of the Hilbert 
space is that of a qubit coupled to a spin N/2, where N 
is the number of TLFs. For the procedure similar to that 
of [ill the oscillation periods with k excitations would be 
given by 2TT/v±^k, where v±^k = \/fc {N + 1 — k) ■ v±. 

V. COLLECTION OF FLUCTUATORS 

We now analyze decoherence of a qubit due to multiple 
TLFs. For this purpose, we introduce an ensemble of 
TLFs with energy splittings e/^„, distributed randomly. 
For each fluctuator n we assume a uniform distribution 
of its energy splitting e/_„ in a wide interval AE, with 
probability density Pn = l/AE. The overall density of 
fluctuators is given by vq = N/ AE, where TV is the total 
number of fluctuators in the interval AE. For simplicity 
we assume all the fluctuators to have the same coupling 
to the qubit v± and the same relaxation rate 7/^. The 
interval is much wider than a single peak, AE 3> v±, 
and the total number of TLFs in the interval is = 
I/O AE > 1. 



We find that the physics is controlled by the dimcn- 
sionless parameter P = vq v±. For D <^ 1 the probability 
for two fluctuators to be in resonance with each other 
is low. Once the qubit is in resonance with one of the 
TLFs, the decay law of the qubit's (ct^) takes the form 
(ID). In this regime we take Fa,, to characterize the de- 
cay. We expect that in most situations the oscillations 
in Eq. ^ will decay fast due to the pure dephasing, 
and one will observe a very fast partial (down to half an 
amplitude) decay of (cr^) followed by further decay with 
rate Tav Thus the relaxation rate is given by a sum of 
many well separated peaks, each contributed by a single 
fluctuator. Since the positions of the peaks are random, 
we expect, for i> ^ 1, a, randomly looking dependence of 
the qubit's relaxation rate on the qubit's energy split- 
ting (and a collection of rare peaks for lower densities, 
1/ <C 1). To characterize the statistical properties, we 
determine in Sec. IV Al the relaxation rate, averaged over 
realizations, and its variance. 

For larger D > 1 the situation changes, as the peaks 
become dense, and the probability to have two or more 
fluctuators in resonance with each other is high. We con- 
clude that it is not reasonable anymore to characterize 
the decay of (cTz) by Tav In this limit the coherent oscil- 
lations turn into much faster relaxation. The excitation 
energy is transferred from the qubit to the TLFs on a 
new, short time scale, ~ {i^v±)~^, which we now call 
the relaxation time. The energy remains in the TLFs for 
much longer time (~ l/7i ) before it is released to the dis- 
sipative baths. Yet, if a strong enough fluctuation of the 
TLFs spectral density occurs, coherent oscillations ap- 
pear again. A set of TLFs almost at resonance with each 
other form an effective strongly coupled fluctuator. The 
decay time of the oscillations is due to the background 
density of TLFs rather than due to the coupling to the 
baths. This could be an alternative explanations for the 
findings of Refs.[l3,[H, [H. In Sec. IV Bl we describe these 
two situations. 



A. Independent fluctuators, i/ <§; 1. 

In this regime the relaxation rate F is given by a sum 
of single-fluctuator contributions. For a given realization 
of the ensemble we obtain 

r(e,) = ^r,,(e, -£/,„) . (28) 

n 

Integrating over the TLF energy splittings e , we obtain 
the average relaxation rate 

(r> = / d^e E ^^(^9 - ^f,n) ~ l( y (29) 

Here d^e = de/^i . . . def^N, and the probability distribu- 
tion p'^) = Yin Pn is given by the product of single-TLF 
distribution functions. p„ = -J^ = 



We also find the variance 



Thus 



((£!)) 1 



(30) 
(31) 



This result can be expected. In the regime v <^ I in 
each realization of the environment the function r(eg) is a 

collection of rare peaks of height ^7/ and width v±. The 
average value of F is, thus, small, but the fluctuations 
are large. As expected, the relative variance decreases as 
the effective density P increases. 

As we have seen in Sec. IIV CI the contribution from 
two TLFs in resonance differs from the sum of two single- 
TLF contributions. This effect leads to modifications of 
Eqs. (p9)) and ([30)) with the further increase in the spec- 
tral density P of the fluctuators. The relaxation rate be- 
comes lower than the result, Eq. (|29p . in the approxima- 
tion of independent fluctuators, and the straightforward 
estimate gives: 



(F) cx 7/ (i^ - cii'^) , 



(32) 



Similarly, for the variance (which is close to the mean 
square) we find 



((F^)) (X (7f)^ {D 



C2V 



(33) 



Here both prefactors ci,C2 ~ 1. For example, ci = 
2-V2,C2 « 1.71 in the rough approximation, when we (i) 
account for correlations by using the rate [Eq. ([26ll ] for 
two resonant TLFs to describe the joint effect of two fluc- 
tuators, n and to, in a certain range around resonance, 
i.e., when their energy splittings differ by less than the 
coupling strength, |e/_„ — €/.,„! < wj^; and (ii) neglect 
correlations for larger detunings |e/.„ — e/,Tn|- 

Figure [8] shows the relaxation rate F for one possible 
realization of the TLF distribution at v = 0.5. The solid 
blue line corresponds to the approximation of indepen- 
dent fluctuators [Eq. (|28|) ]. while the dashed red line is 
calculated using the approximation described. 

We see, that the experimental data, where "random" 
behavior of the relaxation rate as a function of the qubit's 
energy splitting was observed, could be consistent with 
the situation depicted in Fig. [H i.e, with v ^ \. 



B. Spectrally dense fluctuators, v ^ 1 

For higher densities v the calculations above are no 
longer valid. In this section we discuss the limit of very 
high spectral densities, v ^ 1. In the following we re- 
strict ourselves to the one-excitation subspace of the sys- 
tem and neglect the couplings to the baths. Thus we con- 
sider the one-excitation subspace of the following Hamil- 
tonian 

1 1 1 



H 



2 ^ 



(34) 




800 Cg/Tl 



FIG. 8: r as a function of the qubit's level splitting in 
units of the TLF relaxation rate 7/ for one possible realiza- 
tion of the fiuctuator distribution with 1? = 0.5 and v± = 67/. 
Solid blue line: approximation of independent fluctuators 
(|28p . Dashed red line: with account for correlations, see text. 
The corrections are most pronounced in areas, where more 
than one fiuctuator is in resonance (F > ^7/). 



Our purpose is to diagonalize the Hamiltonian in the 
one-excitation subspace and to find the overlap of the ini- 
tial state \i) (qubit excited, all fluctuators in the ground 
state) with the eigenstates |(/)fe), labeled by an index k 
and having the eigenenergies Ek- This allows us to ob- 
tain the time evolution of the initial state: 



t>k\i) e 



(35) 



We begin with the case of a completely uniform spec- 
tral distribution. This case is well known in quantum 
optics as the Wigncr-Wcisskopf theory [l^l . We obtain a 
Lorentzian shape of the overlap function. Figure O shows 
the overlap for an effective density = 1 as a 

function of Ek- We arrive at the probability amplitude 
to find the qubit still excited after time t > 0: 



dEj2\{M^)\^ S{E - Ek)e 



-iEt 



dE 



2- {^^y+E^ 



-iEt 



e-^* , (36) 



where F = ^ D ■ v±. Thus the decay of the initially ex- 
cited qubit in this situation is described by a simple ex- 
ponential decay, | {i\i{t)) [^ = e^^*. The width F of the 
Lorentzian in Fig. [9] determines the decay rate of the ex- 
cited state. 

Note, that we did not include here the coupling of ei- 
ther the qubit or the fluctuators to the dissipative baths. 
These couplings will broaden each of the eigenstates by 
an amount ^ j(. As long as this broadening is smaller 
than the resulting decay rate T = ^ D ■ v± (for strongly 
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FIG. 9: (a) Level structure of a qubit coupled to a uni- 
form spectral distribution of fluctuators. The one-excitation 
subspace is shown. The state, where the qubit is excited, 
|e ttt • • •), is coupled with strength v± to all other levels in 
this subspace. The ground state \g ||| . . .) is energetically 
well separated from the subspace with one excitation, (b) 
Overlap of the initial state \i) = |e m . . .) with the eigen- 
states of the coupled system \(l)k) for a uniform distribution 
of TLFs with an effective density = 1 as depicted above. 
The energy is counted from the energy of the initial state. 



coupled fluctuators {v± ^ 7/) and D > 1 this is al- 
ways the case), the dissipative broadening has little ef- 
fect. Thus the system of the coupled qubit and fluc- 
tuators remains for a long time (~ l/7i) m the one- 
excitation subspace, but the qubit relaxes much faster 
and the energy resides in the fluctuators. 

Further, we analyze the situation with a large num- 
ber of TLFs, whose energy splittings arc accumulated 
near some value. For instance, this behavior may orig- 
inate from the microscopic nature of the fluctuators. 
As we have seen above, a collection of resonant TLFs 
is equivalent, within the single-excitation subspace, to 
one effective TLF with a much stronger coupling to the 
qubit. This results in two energy levels, separated by this 
new strong effective coupling constant. As we discussed 



0.00 h 



- 50 
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FIG. 10: (a) Level structure for a non-homogeneous spec- 
tral distribution of TLFs. At energies near the qubit's level 
splitting eg the density is increased. The state |e |tt . . .) is 
coupled equally to each level in the one-excitation subspace. 
(b) Overlap of the initial state with the eigenstates of the 
coupled system with a higher local density viocai = 10 be- 
tween e = — 10 and e = 10 (in units of v±). Out of resonance 
the effective density is i> = 1. 



above, this may be the origin of the visible properties 
of strongly coupled TLFs. To illustrate this setting, we 
show typical numerical results in Fig. (TUl The data are 
shown for the situation, where instead of many resonant 
TLFs we have a large collection of TLFs distributed in a 
certain energy range. On top of the homogeneous distri- 
bution with density D ~ 1, we assume, locally, a higher 
density {D = 10) of fluctuators with energies close to 
that of the qubit. We obtain a double-peak structure for 
the overlap Performing again a calculation of 

Eq. (j36p . we obtain oscillations with frequency given by 
the energy splitting of the two peaks in Fig. [TOl The 
widths of the peaks (set at this level by the local den- 
sity) determine the decay rate of the oscillations. If these 
peaks are wider than the dissipative broadening, the lat- 
ter can be neglected. Thus, the effect of the fluctuator 
bath with strong density variations on the qubit is cquiv- 
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alent to that of a single TLF with an effective coupHng 
strength v±, much stronger than the couphngs v± n be- 
tween the qubit and the individual physical TLFs. 

VI. CONCLUSIONS 

In this paper we examined the effect of strongly cou- 
pled two-level fluctuators on the dissipative dynamics of 
a qubit. Wc have described the following phenomena: 

(a) If the qubit and TLF arc close to resonance, one 
should observe coherent oscillations between them p^ . 
We have analyzed the effects of dissipation and the re- 
laxation rates. These results may be important in the 
studies of single two-level fluctuator systems, e.g., with 
the idea to use them as quantum memory. 

(b) The situation of a qubit coupled to several TLFs 
with degenerate level splittings is equivalent to the cou- 
pling of the qubit to a single effective TLF with a renor- 
malized, stronger coupling strength. 

(c) Collections of TLFs with spectral density v ^ 1 
could show a seemingly random dependence of the qubit 's 



relaxation rate on its level splitting. 

(d) In the dense hmit, u ^ 1, we conclude that uni- 
form distributions of the TLF energies lead to the expo- 
nential relaxation of the qubit. Local strong fluctuations 
of the density in the energy distribution can, however, 
from the qubit 's viewpoint appear as a single strongly 
coupled effective TLF. We conclude that the signatures 
of strong coupling to single two-level systems often ob- 
served in qubit spectroscopy could well arise from weak 
coupling to many nearly resonant TLFs. We emphasize 
that this equivalence holds only for the initial state, in 
which the qubit is excited and the fluctuators are not. 
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